Self-similarity of single-channel transmission for electron transport in nanowires 
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We demonstrate that the single-channel transmission in the resonance tunneling regime exhibits 
self-similarity as a function of the nanowire length and the energy of incident electrons. The self- 
similarity is used to design the nonlinear transformation of the nanowire length and energy which, 
on the basis of known values of transmission for a certain region on the energy-length plane, yields 
transmissions for other regions on this plane. Test calculations with a one-dimensional tight-binding 
model illustrate the described transformations. Density function theory based transport calculations 
of Na atomic wires confirm the existence of the self-similarity in the transmission. 



I. INTRODUCTION 

When a nanowire is attached to two electron reservoirs 
with different chemical potentials, the electric current 
flows through it. Studies of the electrical conductance 
at the nanoscale have resulted in discoveries of interest- 
ing physical effects, which include conductance quanti- 
zation, molecular rectification, negative differential re- 
sistance, hysteresis, oscillatory length-dependence of the 
conductance^ These and other unusual and potentially 
useful properties make nanowires promising candidates 
for creating various electronic nanodevices^ 

Electrical conduction in nanowires can be approxi- 
mated under certain conditions as a ballistic scattering of 
electrons along the one-dimensional quasiperiodic super- 
lattice of finite length. Since the first fabrication of quasi- 
periodic GaAs — AlAs multilayer structure^ propaga- 
tion of elementary excitations (electrons, photons, polari- 
tons, spin waves) through one-dimensional fractal struc- 
tures has been extensively studied^ It has been demon- 
strated, in particularly, that the quasi-periodic structure 
of the medium manifests itself in the self-similarity of the 
transmissioniSiLSi^ii . The Landauer conductance of the 
generalized Thue-Morse and Fibonacci one-dimensional 
lattices has also been investigatediiiiiSiiiiiiiSiiSiii In ad- 
dition, the fractal conductance fluctuations have been 
predicted and measured in two-dimensional electronic 
billiards in magnetic fields. 18 If the system has a 
quasiperiodic structure and consists of a number of 
identical or similar scattering blocks, the transmission 
through a single block can be mapped into the transmis- 
sion through the entire system. This is the central idea 
behind several powerful approaches to the calculation of 
transmission through quasiperodic structures, those like 
transfer matrix 19 and Dyson-equationsSi The aim of our 
paper is to explore the possibility of the self-similarity 
for single-channel transmission for electron transport in 
nanowires. 

We demonstrate that a single-channel transmission in 
the resonance scattering regime exhibits self-similarity 
as a function of the incident electron energy and the 
nanowire length. We use one-dimensional tight-binding 
model with nearest-neighbor hopping and Landauer the- 
ory to explain the self-similarity of the transmission co- 
efficient for electron transport in nanowires. Based upon 



resonance condition for the transmission, we derive ana- 
lytical formulas for nonlinear continuous transformation 
of the wire length and the energy of the transmitted elec- 
trons. The transformation preserves the number of trans- 
mission maxima by deforming initial region on the energy 
- wire length plane. We confirm the existence of the self- 
similarity by numerical calculation of the transmission for 
various wire lengths, as well as by comparison with the 
first principle transport calculations for sodium wires. 



II. SELF-SIMILARITY OF TRANSMISSION 

To describe a nanowire of length N, we adopt the stan- 
dard tight-binding model with the nearest neighbor con- 
stant hopping rate V and the site energy Eq. The sys- 
tem Hamiltonian is thus described by the three-diagonal 
N x N matrix: 
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The interaction between the wire and the leads is not 
considered explicitly and is accounted for by the intro- 
duction of the self-energy matrix, into the wire 
Green's function. 



G(E) = [i7-H-E(J5)] _1 . 



(2) 



Assuming that the molecule is connected to the left 
(right) lead through its l(iV)th sites, the transmission 
is determined by the IN element of the Green's function 
©: 21 ^ 



T(E) = 4A 2 (E)\G 1N (E)\ 2 , 



(3) 



where A(E) is the imaginary part of the self-energy func- 
tion 



E(J5) = A(E)-iA(E). 



(4) 
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The conductance at low temperature and small voltage 
is given (in units of e 2 /irK) by the Landauer formula 

9 = T(E F ), (5) 

Ep being the Fermi-energy of the electrodes. 

The matrix element of the molecular Green's func- 
tion G\n{E) can be analytically computed for our model 



Hamiltoniani2ii22 To simplify the presentation, we define 
the self-energy within the broadband approximation as- 
suming that A = in eq.Q and considering A as an 
energy- independent constant M> The assumption can be 
relaxed and its consequences are discussed later. 

The analytical calculations yield the following expres- 
sion for the transmission^ 



T(E) = 4£ 2 



sin{0} 

sin{(7V + 1)0} + 2z£sin{iV0} - £ 2 sin{(iV - 1)6*} 
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where we have introduced the dimcnsionlcss quantities 
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If 

eqs. 



E—E, 



tt— I > 1, then the trigonometric function in 
are transformed into the corresponding hyper- 
bolic functions. This results in the exponential scaling 
of the transmission with the nanowire length Nm^ If 
E 2V° I — 1> then eq.© describes the resonance tunnel- 
ing regime where the transmission oscillates as a function 
of the nanowire length A. Oscillatory regime is the in- 
trinsic property of the single channel atomic nanowires 
such as Au and Na; 23 i 24 i 25 i 26 i 27 i 28 to which the one- 
dimensional tight-binding Hamiltonian Q is applicable. 
Therefore we restrict our consideration to this regime and 
demonstrate that the transmission exhibits scaling and 
self-similarity. 

Suppose that the coupling between the wire and the 
lead is weak, i.e. £ < 1. The transmission T(E) ~ £ 2 + 
0(£ 3 ) is very small for £ <C 1 unless sin{(iV+ 1)0} in the 
denominator of eq.© tends to zero. If we put sin{(A" + 
1)0} = 0, which is equivalent to assuming that 

° k= ATI' E k = E o + 2Vcos{0 k }, k = l,2,...,N 

then the transmission achieves its maximum 



T(E k ) 



1 



i+e cos 2 {0 fc }' 



(9) 



For any given length of the wire there are TV different 
resonances © in the transmission. The resonance label 
k in eqs.©, © c an be interpreted as the number of the 
transmission maxima achieved for the given molecular 
length A within the energy interval [Ek , Eq + 2V) . 

The resonance condition © suggests the existence of 
the self-similarity of the transmission T(E) as we vary 
the nanowire length N and energy E of transmitted elec- 
trons. Indeed, the following simple linear transformation 



0/3, J > 1 



fll) 



leaves the resonance condition © unaffected, thereby 
preserving the number of resonances of the transmission 
by deforming the energy interval. The mapping parame- 
ter j = (N + 1)/(N+ 1) is not, in general, an integral but 
a rational number. If the angles and are expressed 
explicitly in terms of the corresponding energies E and 
E via eq. Q , then the transformation i|ll|) becomes non- 
linear: 



E = E + 2V r cos(arccos 



E-E 
2V 



/j), J > 1. (12) 



N=j(N + l)-l, 



(10) 



The transformation (|l(Jf> . l|12fl leads to the self-similarity 
of transmission at [E, E + 2V] and [E, E + 2V] energy 
intervals: the wire of the length A^ within the interval 
[E, E + 2V] has exactly the same number k of the trans- 
mission maxima as the wire of the length A within the 
interval [E,Eq + 2V\. Furthermore, consider the trans- 
mission T(E) of a wire of the length N when the energy of 
incident electrons varies within the range [E±, E%\. Then 
T(E) can be mapped into T(E) for any N > N within the 
window \Ei,E2\, provided that j is defined via eq. l|10fl . 
The number and positions of the conductance peaks © 
for T(E) within [E\, Eq\ will be exactly the same as that 
of f (E) within [E U E 2 ]. 

III. EXAMPLE CALCULATIONS 

A. Numerical calculations with one-dimensional 
tight-binding model 

We confirm the existence of the self-similarity of a 
single-channel transmission by numerical calculations via 
eq.©. The computations are performed for Eq — eV, 
V = 2.5 eV and £ = 0.04. Several illustrative results 
are presented in Fig.l for rational mapping parameters 
j. The transmission T(E) is computed within the energy 
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then the matrix S reduces to the three-diagonal form: 




FIG. 1: Transmission vs. energy for N — 21 (upper panel), 
34 (middle panel), and 107 (bottom panel). For N = 21, 
the energy window 1.5 < E < 4.5 has been chosen. For 
N = 34 (j = 1.59) and 107 (j = 4.91) the corresponding 
energy windows are 3.50 < E < 4.80 and 4.83 < E < 4.98 



windows calculated as prescribed by eas. (|10fl and (|12|l . 

The transmissions for N — 21, 34 and 107 exhibit almost 
perfect self-similarity. Note that the scale of the vertical 
axes in the figures is not the same. This means that 
the transmission curves should be additionally stretched 
in order to make a perfect match between them. Fig. 
1 shows that the widths of the peaks are broader for 
shorter nanowires. This has the following explanation. 
If we expand transmission © in the vicinity of E w Ef. 
and retain the leading on E — Ek and £ contributions, we 
obtain: 



T(E) 



1 



NV(E-E k ) 



(13) 



Thus, the higher is N and/or k, the narrower is the 
Lorentzian peak-shape (|13fl . 

The procedure described above opens the possibility 
of two-dimensional continuation of transmission on the 
E — N plane. The transmission T(E) in the rectangular 
domain [Ei < E < E 2 ], [Ni < N < N 2 ] can be mapped 
into that of the domain [E x < E < E 2 ], [N\ < N < N 2 ] 
via the transformation IjlQI I12|) . 

Most tight-binding and first principle transport calcu- 
lations relay on the nonorthogonal atomic basis functions 
to compute the wire Green's functions. 29 Within the non- 
orthogonal representation the Landauer formula © re- 
mains valid, but the definition of the Green's function © 
should be modified according to 



G(E) = [ES H Y.{E)}-\ 



(14) 



where S is the overlap matrix. If we assume that only 
the nearest-neighbor basis functions overlap significantly, 
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(15) 



where |A| < 1 is the overlap integral between the nearest 
neighbor basis functions. In this case, the expression for 
the transmission © remains correct if the parameter 9 
is redefined: 



cos# = 



E — E a 
2{V-EX)' 



(16) 



Therefore, the non-orthogonality of the basis functions 
breaks the mirror symmetry of T(E) with respect to the 
sign of the site energy Eq, and different choices of the 
site energy Eq are no longer equivalent. 

Having adopted new definition (|16fl for the parameter 
8, we observe that the self-similarity transformation IjlQI 
I11J) remains unchanged. The counterpart of the transfor- 
mation l|12l) in the energy domain should be modified to 
account for the basis nonorthogonality: 



E 



E Q + 2ycos(arccos{ 2( ^° A) } /j) 
1 + 2A cos(arccos { 2 $~_%) } If) ' 



i>i. (17) 



If A = 0, then the transformations i|12|) and i|17fl coincide, 
evidently. If A ^ 0, the non-orthogonality of the basis 
functions causes, through the (energy-dependent) scal- 
ing factor (1 + 2A cos(0/j)) -1 , an additional stretching 
(shrinking) in the energy space, depending on whether 
the product Xcos(0/j) is negative (positive). For exam- 
ple, if we assume that the overlap parameter A = 0.2 
and plot the transmissions analogies to those in Fig. 1, 
we obtain exactly the same curves. However, according 
to eas. <|10|) and (|17f) . the energy intervals will shrink (we 
use eV units for energy E): 1.34 < E < 3.31 for N = 21, 
2.73 < E < 3.47 for N = 34, and 3.49 < E < 3.56 for 
N = 107. On the contrary, if we put A = —0.2 (this is 
equivalent to the consideration of T(—E) for A = 0.2) 
and plot the same transmissions, the corresponding en- 
ergy intervals will stretch as compared with those in Fig. 
1: 1.70 < E < 7.03 for N = 21, 4.86 < E < 7.79 for 
N = 34, and 7.88 < E < 8.27 for N = 107. 



B. First principle simulations for sodium nanowires 

In this section, we examine the prediction of trans- 
mission self-similarity by performing the first princi- 
ple simulations of transport properties of Na atomic 
wires. The calculations were performed using our 
plane wave/pseudopotential implementation of the den- 
sity functional based non-equilibrium Green's function 
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techniques in the CPMD packagei^Mi Our implemen- 
tation is based on nonorthogonal Wannier-type atomic 
orbitals. 30 All systems were treated employing periodic 
boundary conditions and the Kohn-Sham orbitals were 
expanded in plane waves (50 Ry cutoff) at the T point 
of the Brillouin zone. We used local density approxima- 
tion for the exchange-correlation functional and Stumpf, 
Gonze, and Schettler pseudopotentials^ for core elec- 
trons. The system is treated as a "supermolecule" placed 
in a large supercell. The size of supercell is chosen in such 
a way that the distance between the nearest atoms in the 
neighboring cells is larger than 8.5A, so that the inter- 
action between supercell images is negligible. The whole 
system is divided into three parts: left electrode, central 
wire, and right electrode. The electrode part is obtained 
by cutting a few atoms from Na (001) surface. The ge- 
ometry of the electrodes is fixed to the bulk values. The 
wire part is a single chain of Na atoms, where the distance 
between the atoms is constrained to the nearest neighbor 
distance in the bulk system. The distance between the 
electrode part and the wire part is optimized. 



The results of our first principle transport calculations 
are illustrated by Figs. 2. Figs. 2a and 2b depict the 
transmission calculated for Na nanowire with TV = 3 and 
5, respectively. Figs. 2c and 2d show zoomed portions of 
Figs. 2a and 2b, which are seen to exhibit quite a good 
resemblance. 



To link these calculations to the theory developed in 
the paper, the parameters of the tight-binding molecu- 
lar Hamiltonian H, the self-energy matrix 1j(E) and the 
overlap matrix S were extracted from the first principle 
transport calculations. Since we use nonorthogonal ba- 
sis set, the Green's function should be taken in the form 
(I14jl which accounts for the overlap matrix S between 
Wannier type atomic orbitals. The calculated molecu- 
lar Hamiltonian and overlap matrix are described fairly 
well by the three-diagonal matrices and 1(15(1 , respec- 
tively, with the parameters Eq = —2.88 eV, V = 1.25 
eV, and A = —0.28. If we insert these parameters into 
eqs. (110(1 and (|17[1 . then the energy window [—0.56,0] for 
N = 3 will be mapped into [-0.57, -0.35] for N = 5. 
This correlates reasonably well with the actual window 
[—0.64,-0.35]. Of course, the calculated H and S ma- 
trices comply with simple three-diagonal formulas 
and 1|15|) only approximately, and the actual self-energy 
function T,(E) = A(E) — iA(E) is neither purely imagi- 
nary, nor energy independent. Evidently, the presence of 
A(E) causes an additional shift of transmission max- 
ima, and E(-EJ) brings an additional energy dependence of 
T(E) which has not been taken into account within the 
present theoretical analyses. If, however, £(-E) does not 
change significantly within the interval[£i < E < EJ\, 
then the self-similarity in T(E) should be preserved, as 
is proven by Figs. 2c and 2d. 




E,e. V. E,e. V. 



FIG. 2: First principle simulations for Na nanowire. Trans- 
missions calculated for N = 3 and N = 5 are presented in 
Figs, (a) and (b), respectively. Fig. (c) depicts a portion of 
Fig. (a) at the energy interval [—0.56,0], and Fig. (d) shows 
a part of Fig. (b) at the interval [-0.64, -0.35]. The Fermi 
energy is set to zero. 



IV. DISCUSSION AND OUTLOOK 

In this paper, we have demonstrated that a single- 
channel transmission for electron transport in nanowires 
exhibits the self-similarity in resonance tunneling regime. 
Starting from the resonance condition, we proposed the 
transformation of the wire length and the incident elec- 
tron energy, which preserves the number of maxima 
in the transmission. A sample calculations on one- 
dimensional tight binding model produces exact agree- 
ment with analytical results. The density functional 
based transport calculations of Na atomic wires confirm 
the existence of self-similarity in the transmission for re- 
alistic system. 

The existence of scaling properties and self-similarity 
of the molecular transmission is also supported by a pos- 
sibility to define it via the Friedel sum rule 33 . It is not 
necessary to have a detailed information on a large va- 
riety of the parameters describing the molecule and the 
leads to predict transport properties of a single channel 
conductor. To compute the single-channel conductance 
by the Friedel sum rule, it is sufficient to know only the 
difference between the number of electrons occupying the 
even and the odd eigenstates of the molecule. 

It is thus not very unrealistic to suggest that the 
self-similarity in molecular transmission and conductance 
may be a rather general phenomenon. In terms of pos- 
sible applications, our results suggest that the transmis- 
sion, computed within a certain energy and length win- 
dow, can be transferred to other windows via the trans- 
formation 1(10(1. I|12|l . I|17|) . For example, the transmission 
of a short wire within a relatively wide energy window 
can be mapped into the transmission of a longer wire 
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within a narrower energy window. If we assume, ad- 
ditionally, that the site energy E can be purposefully 
modified )? 4 ^ then the molecular conductance g, eq.JHJ), 
becomes a function of the (modified) site energy Eq and 
molecular length N and also exhibits the self-similarity. 
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